Resonant-state solution of the Faddeev-Merkuriev integral equations for three-body 

systems with Coulomb potentials 
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A novel method for calculating resonances in three-body Coulombic systems is proposed. The 
Faddeev-Merkuriev integral equations are solved by applying the Coulomb-Sturmian separable ex- 
pansion method. The e~e^e~ S-state resonances up to n = 5 threshold are calculated. 



I. INTRODUCTION 

For three-body systems the Faddeev equations are the 
fundamental equations. Three-body bound states cor- 
respond to the solutions of the homogeneous Faddeev 
equations at real energies, and resonances, as is usual in 
quantum mechanics, are related to complex-energy solu- 
tions. 

The Faddeev equations were derived for short-range 
interactions. However, if we simply plug-in a Coulomb- 
like potential they become singular. A formally exact 
approach was proposed by Noble 0. His formulation 
was designed for solving the nuclear three-body Coulomb 
problem, where all Coulomb interactions were repul- 
sive. The interactions were split into short-range and 
long-range Coulomb-like parts and the long-range parts 
were formally included in the "free" Green's operator. 
Merkuriev extended the idea of Noble by performing the 
splitting in the three-body configuration space This 
was a crucial development since it made possible to treat 
attractive Coulomb interactions on an equal footing with 
repulsive ones. 

Recently we have presented a method for treating 
the three-body Coulomb scattering problem by solv- 
ing Faddeev-Merkuriev integral equations using the 
Coulomb-Sturmian separable expansion technique 
We solved the inhomogeneous Faddeev-Merkuriev inte- 
gral equations for real energies. Previously, for calculat- 
ing resonances in three-body systems with short-range 
plus repulsive Coulomb interactions, we solved homo- 
geneous Faddeev-Noble integral equations by using the 
Coulomb-Sturmian separable expansion technique Q . In 
this paper by combining the concepts of Refs. jsj and 
we solve the homogeneous Faddeev-Merkuriev integral 
equations for complex energies. This way we can han- 
dle all kind of Coulomb-like potentials in resonant-state 
calculations, not only repulsive but also attractive ones. 

In section II we present the homogeneous Faddeev- 
Merkuriev integral equations, outlined for systems where 
two particles out of the three are identical. Many sys- 
tems, like e~e"*"e~ and H~ , fall into this category. Then, 
in section HI, we present the solution method adapted to 
the case where all charges have the same absolute value. 
In section IV we present our calculations for the L — 



resonances of the e^e^e^ system up to the n — 5 thresh- 
old and compare them with the results of complex scaling 
calculations Q. 



II. FADDEEV-MERKURIEV INTEGRAL 
EQUATIONS 

The Hamiltonian of a three-body Coulombic system 
reads 



H = H 







(1) 



where i?" is the three-body kinetic energy operator and 
denotes the Coulomb-like interaction in the subsys- 



tem a. We use throughout the usual configuration-space 
Jacobi coordinates Xa and i/a. Thus v, 
Xa — Va{xa)). The Hamiltouian 
three-body Hilbert space. The two-body potential oper- 
ators are formally embedded in the three-body Hilbert 
space 



'-^ only depends on 
1) is defined in the 



{x)ly, 



(2) 



where \y is a unit operator in the two-body Hilbert space 
associated with the y coordinate. We also use the nota- 
tion X = {XaTUa} G R^. 

The role of Coulomb potentials in Hamiltonian (|l|) are 
twofold. Their long-distance parts modify the asymptotic 
motion, while their short-range parts strongly correlate 
the two-body subsystems. Merkuriev introduced a sepa- 
ration of the three-body configuration space into different 
asymptotic regions. The two-body asymptotic region 
is defined as a part of the three-body configuration space 
where the conditions 



\x^\<xl{l + \v^\/vlf/'^, 



(3) 



with xi^,y^ > and v > 2, are satisfied. Merkuriev 
proposed to split the Coulomb interaction in the three- 
body configuration space into short-range and long-range 
terms 



,(0 



(4) 



1 



where the superscripts s and I indicates the short- and 
long-range attributes, respectively. The splitting is car- 
ried out with the help of a splitting function (a which 
possesses the property 



1, Xc^en a 

otherwise. 



(5) 



In practice, in the configuration-space differential equa- 
tion approaches, usually the functional form 

C(x, y) = 2/(1-1- exp [{x/x"r/{l + y/y")] } , (6) 

was used. 

The long-range Hamiltonian is defined as 



,(0 ^,,(0 



"2 ^ ''3 ' 

and its resolvent operator is 

GW(z) = (z-i?W)-\ 



(7) 



(8) 



where z is the complex energy-parameter. Then, the 
three-body Hamiltonian takes the form 



H = i/W 



(9) 



which formally looks like a three-body Hamiltonian with 
short-range potentials. Therefore the Faddeev method is 
applicable. 

In the Faddeev procedure we split the wave function 
into three components 



I*) = IV-l) + IV'2) + IV'3), 

where the components are defined by 



(10) 



(11) 



In case of bound and resonant states the wave-function 
components satisfy the homogeneous Faddeev-Merkuriev 
integral equations 



\i;,)=Gf\z)v['\\i;2) + m] 
\^2)=G^^\z)vi'\\i^,) + \^,)] 

|V'3)=4'^(z)4^^[|^i) + |^2)] 



(12) 
(13) 
(14) 

at real and complex energies, respectively. Here Gq'' is 
the resolvent of the channel long-ranged Hamiltonian 

^ H^'^ + vi'\ (15) 



Merkuriev has proved that Eqs. 



(p^HlJ) possess compact kernels, and this property re- 
mains valid also for complex energies z — E~iT /2, F > 0. 

In atomic three-particle systems the sign of the charge 
of two particles are always identical. Let us denote them 
by 1 and 2, and the non-identical one by 3. In this case 



is a repulsive Coulomb potential which does not sup- 
port two-body bound states. Therefore the entire v'^ can 
be considered as long-range potential. The long-range 
Hamiltonian is modified as 



„(0 



,(0 



"3 • 



Then, the three-body Hamiltonian takes the form 



H = ifW 



is) , is) 



(16) 



(17) 



i.e. the Hamiltonian of the system appears formally as a 
three-body Hamiltonian with two short-range potentials. 
Therefore the Faddeev procedure, in this case, gives a set 
of two-component Faddeev-Merkuriev integral equations 



IV^i) 



1^2) = G 



IV-i). 



(18) 
(19) 



Further simplification can be achieved if the particles 1 
and 2 are identical. Then, the Faddeev components 
and IV'2)) in their own natural Jacobi coordinates, have 
the same functional form 



(a^iyilV'i) = (a;22/2|V'2) = {xy\i!). 



(20) 



Therefore we can determine IV') from the first equation 
only 



IV') 



(21) 



where V is the operator for the permutation of indexes 
1 and 2 and p — ±.1 are eigenvalues of V. We note that 
although this integral equation has only one component 
yet gives full account on asymptotic and symmetry prop- 
erties of the system. 



III. SOLUTION METHOD 

We solve these integral equations by using the 
Coulomb-Sturmian separable expansion approach 
The Coulomb-Sturmian (CS) functions are defined by 



{r\nl) 



(n + 2; + l)! 



1/2 



(26r)'+i cxp(-6r)L2'+i(26r), 

(22) 



with n and / being the radial and orbital angular mo- 
mentum quantum numbers, respectively, and b is the 
size parameter of the basis. The CS functions {|?^/)} 
form a biorthonormal discrete basis in the radial two- 
body Hilbert space; the biorthogonal partner defined by 
{r\nl) = {r\nl)/r. Since the three-body Hilbert space is a 
direct product of two-body Hilbert spaces an appropriate 
basis can be defined as the angular momentum coupled 
direct product of the two-body bases 

\nvlX)a = \nl)a(it)\vX)a, (n, = 0, 1, 2, . . . ), (23) 
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where \nl)a and \i^X)a are associated with the coordinates 
Xa and Ha, respectively. With this basis the complete- 
ness relation takes the form (with angular momentum 
summation implicitly included) 



lim 



N 
riM—O 



ni'lX)c 



(nvlW 



lim 1 

N^oo 



N 



(24) 



Note that in the three-body Hilbert space, three equiv- 
alent bases belonging to fragmentation 1, 2 and 3 are 
possible. 

We make the following approximation on the set of 
Faddeev-Merkuriev integral equations 



|^l) = Gf'(z)lf.i^)[lf|^2) + l^|V'3)] 
|^2)=G«(z)l^(^)[lf|^l) + l^|V'3)] 
|^3)=G^"(^)lf4^'[lf|V^l) + l^|V'2)], 



(25) 
(26) 
(27) 



i.e. the short-range potential Va in the three-body 
Hilbert space is taken to have a separable form, viz. 



N 



^0 — 



J2 \n^l^)c.v!S pin'i^'l'X'l, (28) 



where G^^ = i{nh'lX\Gi'\n'iy'l'X')i. A unique solution 
exists if and only if 



.(0 



det{[G(')(z)] 



m } 



(32) 



Unfortunately Gi'' is not known. It is related to the 
Hamiltonian h[''\ which itself is a complicated three- 
body Coulomb Hamiltonian. In the three-potential for- 
malism 1^ Gi'^ is linked to simpler quantities via solution 
of a Lippmann-Schwinger equation. 



{Gf^r'^{G,r'-u„ 



where 



and 



{nvlX\Gi\n'iy'l'X')i 



C^i„./A,„'.'rA' = i(«^^A|C/i|nVrA')i. 



(33) 



(34) 



(35) 



In our special case, where the sum of the charges of par- 
ticles 2 and 3 is zero, the operator Gi is the resolvent 
operator of the Hamiltonian 



n,u,n' ,i''—0 



Hi = H" 



(36) 



where v^^l = a{wlX\vi'^\n'i^'r X')p. In Eq. (H) the ket 



and bra states are defined for different fragmentation, 
depending on the environment of the potential operators 
in the equations. The validity of this approximation re- 
lies on the square integrable property of the terms like 
Va ^\ipf3)i which is guaranteed due to the short range na- 
ture of Va \ 

For solving Eq. ^W) we proceed in a similar way. 



_ /^(OnAr„(s) 



G\'>r^vrpPl^\i;), 



(29) 



(s) 

i.e. the operator v{ pV in the three-body Hilbert space 
is approximated by a separable form, viz. 



lim V^vY'pVl\ 

N~>oo 

l^vi'^pVl^ 

N 



^ \nvlX)iyXi{n'v'l'X' 



(30) 



n,u.n'u'—Q 



where v_i 



= i{nvlX\v['\v\n'v'l' X')i. Utilizing the 
properties of the exchange operator V these matrix 
elements can be written in the form -ir^ = p x 

i{niylX\v['^\n'p'l'y)2. 

With this approximation, the solution of Eq. ( |2l| ) turns 
into solution of matrix equations for the component vec- 
tor ip = i{nh'lX\ipi) 



}t, = 0, 



(31) 



and the polarization potential Ui is given by 

U,=v!i^+v§. (37) 

The most crucial point in this procedure is the cal- 
culation of the matrix elements Gi, since the potential 

matrix elements Vi"^ and lJ_i can always be evaluated nu- 
merically by making use of the transformation of Jacobi 
coordinates 0. The Green's operator Gq is a resolvent of 

hy 

C(^.\ „„A h — f^o^^ which act in 
different two-body Hilbert spaces. Thus, according to 
the convolution theorem the three-body Green's oper- 
ator Gq, equates to a convolution integral of two-body 
Green's operators, i.e. 



the sum of two commuting Hamiltonians, Hi — h^i Tn^yi , 
where h^^ = h'^^ + {xi) and hy-^ 



1 



Gi{z) = —.f dz' g^,{z - z') gy,{z'), 



(38) 



where 5^1(2) = [z - h^i) ^ and gy^iz) = (z - hyj ^ 
The contour C should be taken counterclockwise around 
the continuous spectrum of hy-^ such a way that gx^ is 
analytic on the domain encircled by C. 

To examine the structure of the integrand let us shift 
the spectrum of gx^ by taking z = E + ie with positive e. 
By doing so, the two spectra become well separated and 
the spectrum of gy-^ can be encircled. Next the contour 
G is deformed analytically in such a way that the upper 
part descends to the unphysical Riemann sheet of gy-^ , 
while the lower part of G can be detoured away from the 
cut [see Fig. 0]. The contour still encircles the branch 



3 



cut singularity of gy-^, but in the e ^ limit it now 
avoids the singularities of g^-^. Moreover, by continuing 
to negative values of e, in order that we can calculate 
resonances, the branch cut and pole singularities of gx^ 
move onto the second Riemann sheet of gy-^ and, at the 
same time, the branch cut of gy-^ moves onto the second 
Riemann sheet of ■ Thus, the mathematical condi- 
tions for the contour integral representation of Gi{z) in 
Eq. can be fulfilled also for complex energies with 
negative imaginary part. In this respect there is only a 
gradual difference between the bound- and resonant-state 
calculations. Now, the matrix elements can be cast 
in the form 

G,{z)^^J dz'gjz-z')gjz'), (39) 

where the corresponding CS matrix elements of the two- 
body Green's operators in the integrand are known an- 
alytically for all complex energies (see |^ and references 
therein), and thus the convolution integral can be per- 
formed also in practice. 

IV. RESONANT STATES IN POSITRONIUM 
IONS 

We calculate resonant states in positronium ion with 
L = total angular momentum. The positronium ion, 
Ps^ or e^e+e", is a three-body Coulomb system that 
consists of two electrons and one positron. We calculate 
its resonances by solving Eq. (pT|). We took a;° = 18ao, 
j/° = 50ao and = 2.1 as the parameters of the splitting 
function, respectively. 

Before presenting our final results we demonstrate the 
convergence properties of this method. In Table (^ we 
show the convergence of a resonant state energy with 
respect to angular momentum channels and number of 
Coulomb-Sturmian basis states N in the expansion. This 
table shows the accuracy and stability of our calculations. 
Table (|ll|) contains the final results. For the low- lying 
resonances we used CS parameter b — 0.25aQ ^, and for 
the high-lying states we took b — O.lSa;^^. We compare 
our calculation with the result of complex scaling calcu- 
lations Ref. §. We can report perfect agreements for the 
position of the resonances, but, in most of the cases, we 
got much smaller values for the width. 

V. CONCLUSIONS 

In this article we have presented a new method for 
calculating resonances in three-body Coulombic systems. 
Our approach is based on the solution of the homoge- 
neous Faddeev-Merkuriev integral equations for complex 
energies. For this, being an integral equation approach, 



no boundary conditions are needed. We solve the in- 
tegral equations by using the Coulomb-Sturmian sepa- 
rable expansion technique. The method works equally 
well for three-body systems with repulsive and attractive 
Coulomb interactions. 
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FIG. 1. Potential v^"' , the short-range part of a —1/x at- 
tractive Coulomb potential. 
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FIG. 2. Potential v^'\ the Ion g- range part of a —1/x at- 
tractive Coulomb potential. 




FIG. 3. Analytic structure of gxi{z — z') gy^{z') as a func- 
tion of z' with z = iJ -I- ie, _E > 0, £ > 0. The contour C 
encircles the continuous spectrum of /ij,^. A part of it, which 
goes on the unphysical Riemann-sheet of Qy^ , is drawn by 
broken line. 
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TABLE I. Convergence of '^S^ 3s4s (L = 0) positronium 
resonance state, b=0.25. 
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TABLE II. Doubly excited L = resonances of Ps . The 
energies and widths are expressed in Rydbergs. 
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